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The discrete spectrum of fluctuations of the metric about an AdS^ black hole background 
are found. These modes are the strong coupling limit of so called glueball states in a dual 
3-d Yang-Mills theory with quantum numbers J^'-' = 2+"'", 1 ^, 0"*"^. For the ground state 
modes, we find the mass relation: m{0^~^) < m{2~^~^) < m(l Contrary to expectation, 
the mass of our new O"*""*" state (m^ = 5.4573) associated with the graviton is smaller than 
the mass of the O"*""*" state (m^ = 11.588) from the dilaton. In fact the dilatonic excitations 
are exactly degenerate with our tensor 2+"'" states. We find that variational methods gives 
remarkably accurate mass estimates for all three low-lying levels while a WKB treatment 
describes the higher modes well. 
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1. Introduction 

The Maldacena conjecture [Q and its furtlier extensions p|0 allow us to compute 
quantities in a strongly coupled gauge theory from its dual gravity description. In particu- 
lar, Witten 1^ has pointed out if we compactify the 4-dimensional conformal super Yang 
Mills (SYM) to 3 dimensions using anti-periodic boundary conditions on the fermions, then 
we break supersymmetry and conformal invariance and obtain a theory that has interest- 
ing mass scales. In Refs. and 0, this approach was used to calculate a discrete mass 
spectrum for O"*""*" states associated with Tr[F^] at strong coupling by solving the dilaton's 
wave equation in the corresponding gravity description. Although the theory at strong 
coupling is really not pure Yang-Mills, since it has additional fields, some rough agreement 
was claimed between the pattern of glueball masses thus obtained and the masses com- 
puted for ordinary 3-d Yang-Mills by lattice calculations. In Ref. the calculation was 
extended to solving the modes of the vector and two form fields with certain polarizations. 
Several other investigations have expanded our knowledge of the relation between string 
theory and gauge theories, in these and related directions . 

Here we consider the discrete spectrum for the graviton field which is dual to the energy 
momentum tensor of the Yang-Mills field theory. This is a very interesting case, since a 5-d 
graviton field has five independent on-shell states that in a dual correspondence to a 4-d 
Yang-Mills theory would account for the five helicity components of a massive 2++ state 
lying on the leading diffractive (Pomeron) trajectory 00. In the present analysis, these 
five states are mapped by dimensional reduction to a reducible multiplet with quantum 
numbers 2+"'", 1 '", 0+"'" in a 3-d Yang Mills theory |]T^. Contrary to expectations, we 
find that the mass gap defined by the lowest 0"'"+ state in the graviton multiplet is lower 
than mass gap for the dilaton. In fact our 2+"'" states are exactly degenerate with the 
lowest dilaton states, which shall be denoted as 0"'""'", while the 1 ^ equation gives the 
most massive states. Remarkably these mass relations, 

m(0++) < m(0++) = m(2++) < m(l-+), 

for the lowest modes are reminiscent of weak coupling glueball spectra from lattice gauge 
theory or the MIT bag model. Extensions to the graviton multiplet for states dual to the 
4-d Yang Mills theory will be reported separately. 

Quite apart from the issue of how closely the gauge theory obtained from the AdS 
description at strong coupling agrees "accidentally" with ordinary Yang-Mills, it is a very 
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interesting fact that one obtains discrete masses representing long lived states in the theory 
from a very different looking calculation done with the gravity variables. On the gravity 
side, if we compactify the AdS^ metric on a large then the metric is a periodic version 
of the original AdS^ x metric. At high temperature (small radius of the S^) this metric 
is replaced by an AdS black hole metric. This geometry has stationary modes for the fields 
of string theory (and in particular the fields of supergravity), and it is conjectured that the 
spectrum of these modes gives the mass spectrum for the compactified gauge theory on 
X S^. Understanding the properties of such a Yang-Mills theory even at strong coupling 
is interesting. 

In this paper we carry out the calculation of the discrete modes for the perturbations 
of the gravitational metric. We restrict the perturbations to be independent of the 
factor in the metric and constant in the co-ordinate for the compact S^. Our results are 
presented as follows. 

In sec II, we derive three different wave equations that arise from fluctuations in the 
gravitational metric. These correspond to spin-2, spin-1 and spin-0 perturbations from the 
viewpoint of the rotation group in the 3-dimensional noncompact space: 

(i) The spin-2 equation can be mapped by a change of variables to the equation for 
a massless scalar, and so the energy levels found are the same as those for the dilaton field 

(ii) The spin-1 equation gives energy levels that are somewhat higher than the spin-2 

case. 

(iii) We first obtain the spin-0 equation in a gauge which has only diagonal compo- 
nents of the metric. We obtain a third order equation for the perturbation. We observe 
that there is a one parameter residual gauge freedom left in the ansatz, and we use the 
form of this gauge mode to reduce the third order equation to a second order equation 
for the perturbation. The situation here is similar to the analysis of polar perturbations 
of the 3-1-1 Schwarzschild hole in a diagonal gauge, where a third order system is also 
found. We comment on the gauge freedom present in the latter case; this discussion may 
illuminate some of the mysteries associated with the existence of the second order 'Zerilli 



equation' [|TT| for the above mentioned polar perturbations of the Schwarzschild metric. As 
a check on our result, we then study the perturbation in an alternative gauge, similar to 
one used by Regge and Wheeler [|12|] . In this latter gauge we directly obtain a second order 



equation, which is seen to be equivalent to the equation obtained by the first method. 



In Sec III, we solve these equations numerically. We comment of the issue of the 
boundary condition to be used at the horizon of the black hole: The physics we impose 
is that there be no flux transfer across the horizon of the black hole, in agreement with 
the boundary condition used in Ref. . (In Appendix B, we also show that the desired 
boundary conditions at r = 1 and r = oo together allow us to treat these equations in a 
Sturm-Liouville approach.) We present accurate numerical solutions to the flrst 10 levels. 
The spin-0 equation gives energy levels that are lower than the spin-1 and spin-2 equations. 

In Sec IV, we consider approximation methods to gain insight into the mass relations. 
The ground state mass is also given in terms of a very simple variational ansatz and the 
entire spectrum approximated very well by a low-order WKB approximation. We also 
demonstrate how the boundary condition at the horizon can be implemented more effec- 
tively with the help of an "effective angular momentum" . (Technical details are postponed 
to Appendix B on the variational method. Appendix C on the Schrodinger form of the 
wave equations and Appendix D on the derivation of the WKB expansion.) 

In sec V, we discuss the results and make a few cautionary remarks on the difficulties 
of the comparison of this strong coupling spectrum with glueballs as computed in lattice 
gauge theory or classified in bag models. We recommend a more thorough analysis of 
the complete set of spin-parity states for the entire bosonic supergravity multiplet and its 
extension to 4-d Yang Mills models. 

1.1. AdS/CFT correspondence at Finite Temperature 

Let us review briefiy the proposal for getting a 3-d Yang-Mills theory dual to super- 
gravity. One begins by considering Type IIB supergravity in Euclidean 10-dimensional 
spacetime with the topology M5 x . The Maldacena conjecture asserts that IIB su- 
perstring theory on AdS^ x is dual to the A/" = 4 SYM conformal field theory on the 
boundary of the AdS space. The metric of this spacetime is 

dr"^ 

ds'^/Rl^^ = r'^idr'^ + dx\ + dxj + dxj) + + dQj , (1.1) 

where the radius of the AdS spacetime is given through -R^^^ = QsNa'"^ {qs is the string 
coupling and Ig is the string length, l"^ = a'). The Euclidean time is r = ixQ. To break 
conformal invariance, following Q , we place the system at a nonzero temperature described 
by a periodic Euclidean time r = T + 2nRo. The metric correspondingly changes, for small 
enough Rq, to the non-extremal black hole metric in AdS space. For large black hole 
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temperatures, the stable phase of the metric corresponds to a black hole with radius large 
compared to the AdS curvature scale. To see the physics of discrete modes we may take 
the limit of going close to the horizon of this large hole, whereupon the metric reduces to 
that of the black 3-brane. This metric is (we scale out all dimensionful quantities) 

ds^ = f{r)dT^ + f-\r)dr^ + r^{dxl + dxj + dxj) + dnj , (1.2) 

where 

= r^-^ (1-3) 

On the gauge theory side, we have a N = A susy theory corresponding to the AdS 
spacetime, but with the compactification we have this theory on a circle with antiperi- 
odic boundary conditions for the fermions. Thus supersymmetry is broken and massless 
scalars arc expected to acquire quantum corrections. 

From the view point of a 3-d theory, the compactification radius acts as an UV cut- 
off. Before the compactification the 4-d theory was conformal, and was characterized by 
a dimensionless effective coupling (gyM)^^ dsN. After the compactification the theory 
is not conformal, and the radius of the compact circle provides a length scale. Let this 
radius be R. Then a naive dimensional reduction from 4-d Yang-Mills to 3-d Yang-Mills, 
would give an effective coupling in the 3-d theory equal to ((7y]^)^A^ = {gYM)'^N/R. 
This 3-d coupling has the units of mass. If the dimensionless coupling (5'y|//)^-^ is much 
less than unity, then the length scale associated to this mass is larger than the radius of 
compactification, and we may expect the 3-d theory to be a dimensionally reduced version 
of the 4-d theory. Unfortunately the dual supergravity description applies at (gyM)^^ » 
1, so that the higher Kaluza- Klein modes of the compactification have lower energy than 
the mass scale set by the 3-d coupling. Thus we do not really have a 3-d gauge theory 
with a finite number of additional fields. 

One may nevertheless expect that some general properties of the dimensionally re- 
duced theory might survive the change between small and large coupling. For this purpose 
we look at the pattern of masses and spins that are obtained for the fields that are singlets 
under the - i.e. we ignore the Kaluza-Klein modes of the S^. In keeping with earlier 
work, we also restrict ourselves to modes that are singlets of the S0{6), since non-singlets 
under the and the S0{6) can have no counterparts in a dimensionally reduced QCD^. 
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2. Wave Equations 

2.1. Field equations and the equilibrium configuration 

We consider Type IIB supergravity in Euclidean 10- dimensional spacetime. The space- 
time will have the topological form M5 x 5"^. For the perturbations that we are interested 
in, we can set to zero the dilaton ^, the axion C and the 2-form fields B^^^^ , B^^. The 
metric and the 4-form gauge field satisfy the Einstein equation 

p.. — 1/? ppark ^2 1^ 

^^flV g flpiJTK-'- O \^'^) 

and the self-duality condition for the field strength of the four form 



TP — —f PMi^P'^-r fool 

^ jiupcTT gi c^jypur/i' i^'p'cr'r' -'^ V / 

Notations: Indices with a 'hat' are 10-d indices. The coordinates of the will be called 
y", and the indices here will be a, /3, . . .. The coordinates of M will be x^, and the indices 
will be /U, z^, . . .. The e tensor is real in spacetime with Lorentzian signature, and thus pure 
imaginary in Euclidean signature. Thus F will be real on 5"^, and imaginary on M5. 
The equilibrium configuration about which we perturb will be given by the metric 

ds^ = f{r)dT^ + f-^{r)dr'^ + r'^{dx\ -f dxl -f dxl) + dVLl (2.3) 

where 

/(r) = r^-^ (2.4) 

The last term in ( |2.3| ) gives the metric of 5"^ and we have chosen units to make this a 
sphere of unit radius. The 5-form field strength with indices in the is a constant times 
the volume form on 5"^. By the self duality condition (^.2|) we have 



^MP-.^r"" = -2A^7m. (2.5) 

where 



2A = ^-F^p^SeF''!'^'' = —Fxp.r.F^P''^'' (2.6) 



is a positive constant. 
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2.2. Ansatz for the perturbations 

We wish to consider fluctuations of tlie metric of the form 



gfiiy = Qiiu + hf^^{x) (2.7) 

Thus the perturbations will have no dependence on the coordinates of the 5"^. Further, 
we wish to keep unchanged all the other fields in the theory. It is easy to see that 
(/), C, B^^^^ J j^RR ^QT^ ]-,g fixed, since they arise quadratically in the action, and 
vanish in the equilibrium configuration. (The metric above is the Einstein metric, and 
so there is no linear coupling of h to the dilaton.) It is also consistent to hold fixed the 
geometry of the and the 5-form field-strength on the S^. Keeping these quantities fixed 
keeps unchanged the flux of F on ^ which is one constraint that the perturbations must 
respect. The self duality condition (2.2) can then regarded as determining the value of F 



on M. But then we see that we will still get after such a perturbation: 
with A the same constant as in ( p.6| ). The Einstein equations then give 



R^, = -^A^^. (2.9) 



Here the Ricci tensor R^^j can be computed solely from the 5-d metric of M; all contri- 
butions involving the variables are zero. Thus if we find a solution to (|2.9| ), we have a 
consistent solution the to entire supergravity problem. With the expansion ( |2.7D we get 

V;A^ - \h\,^.u + \h^,x■,u^ + \Kx,^^ + {D - 1) ^ = (2.10) 

where D = 5. 

To see how many kinds of modes we expect, consider first the graviton in 5-d fiat 
spacetime with coordinates a;^ z = 1, . . . 5, and signature ( — h + + +). Then the graviton 
modes can be brought by gauge freedom to a form 

V=eM-e^'^"'""'^ ^,z^ = 3,4,5, = (2.11) 

Thus there are 5 independent polarizations of the graviton. 

In the spacetime ( |2.3| ) the coordinates xi,X2,xs project to the 3-d spacetime of the 
gauge theory that we wish to study, r is direction of the gauge theory that is compactified in 
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obtaining the 3-d theory from the 4-d theory, and r is the radial coordinate of the spacetime 
that loosely speaking governs the scale at which the gauge theory is being studied. Our 
perturbations will have the form 

V = e/.-We-"'"' (2.12) 

where we have chosen to use xs as a Euclidean time direction to define the glueball masses 
of the 3-d gauge theory. We start with fixing for the gravitational perturbations the gauge 

hs^, = (2.13) 

though later we will have occasion to consider changes to another gauge as well. 

2.3. Wave equations for the perturbations 

From the above ansatz and the metric ( |2.3| ) we see that we have an 5*0(2) rotational 
symmetry in the xi — X2 space, and we can classify our perturbations with respect to this 
symmetry. We then find the following categories of perturbations: 

(a) Spin-2: There are two linearly independent perturbations which form the spin-2 rep- 
resentation of the above SO {2): 

hi2 = h2i = qT{r)e~^^^ , all other components zero (2-14) 

hii — —h22 = qT{r)e~^^'^ , all other components zero (2-15) 

The Einstein equations ( |2.10| ) give, 

1 3 m'^ 1 

^'(1 - ^Wrir) + r(l + ^Jq'Ar) + [— - 4(1 + ^)]?T(r) = 0. (2.16) 

Defining 

^T = % (2.17) 

Eq. ( p.l6|) becomes 

r'il - ^)<j^Ur) + r(5 - ^)Mr) + "^Mr) = <pT;X^ = 0, (2.18) 

which is the free wave equation, and thus the same equation as that satisfied by the 
dilaton (with constant value on the S^). These dilaton modes have been studied 
in 00 for example , and we therefore will obtain the same energy levels for the 
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spin-2 perturbations as those obtained for the dilaton, e.g.^ for the ground state, 
m^ o ~ 11.588. 
(b) Spin-1: 

hir = hri = qv {r)e~'^^^ , i = 1 or i = 2, all other components zero (2.19) 
The Einstein equations (]2.1U[) give, 



1 1 1 

r^{l - :^)q'{.{r) + r(l - :^)q'v{r) + i-;;^ - 4(1 - :;:^)]qv{r) = (2.20) 



(c) Spin-0: Based on the symmetries we choose an ansatz where the nonzero components 
of the perturbation are 

^ii = ^22=gi(r)e-"^-^ 

Kr = -2Qi(r)^e-'^^3 + g2(r)e-"^^3 (2.21) 
Kr = q3ir)e-^^^ 

where /(r) is given in ( |2.4| ). We have chosen to express h^T in such a way so that the 
part involving qi is traceless, which we expect is to dominate at large r. 
The analysis of the Einstein equations are given in Appendix A, where they are derived 
in the above gauge and also in a different gauge which is less directly related to the 



flat space form, (|2.11| ), but more easy to solve: 



/irr = ?i(r)e-"^^^ 

hrr = Q2(r)e-'^^3 (2.22) 
hr3 = hsr = ?3(r)e-"^^^ 



The field equation is, for ^3 = q's(r) 



where 



d 

{^^2(0^ +pi(r)— +po(r)}Q5(r) = 0, (2.23) 



P2(r) = r^{r^ - l)^[3(r^ - 1) + mV^] 

pi(r) = r(r^ - l)[3(r^ - l)(5r^ + 3) +m^r'^{7r^ + 5)] (2.24) 
po(r) = 9{r^ - 1)3 + 2mV2(3 + 2r^ + 3r^) + mV(r^ - 1). 
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3. Numerical Solution 



To calculate the discrete spectrum for our three equation, one must apply the correct 
boundary conditions at r = 1 and r = oo. /,From the indicial equations, the asymptotic 
value of the two linearly independent solutions at infinity are 

qv{r) ~ r"^ (3.1) 
qsir) ~ r"^, r"^ , 

and at r = 1 are 

qrir) ~ 1, log{r - 1) , 

gv(r) ~r-l, {r - l)logir - 1) , (3.2) 
qs{r) {r - 1)"\ (r - l)~^log{r - 1) . 

In all cases the appropriate boundary condition at r = 1 is the one without the 
logarithmic singularity. At r = oo the least singular boundary is required to have a 
normalizable eigenstate. Matching these two boundary conditions results in a discrete set 
of eigenvalues m^, where n is the number of zeros in the wave function inside the interval 
r G (l,cx3). We solved the eigenvalue equations by the shooting method, integrating 
from ri ^ 1 to large Too — oo. The parameter is adjusted by applying Newton's 
method to satisfy the condition lim^^oo q{f) = 0. To begin the integration one needs the 
function and its first derivative at r = 1. Thus expanding to the next significant order 
for the regular solution at r = 1, the initial conditions for the tensor were taken to be 
qriri) = 1, qxiri) = 2 — m? /A and for the vector to be qv{fi) = 0, q'y{ri) = 1. For the 
scalar equation, we found it convenient to define an new function, /s(r) = (r^ — l)qs{r)/r, 
so that the initial value is nonsingular and only the correct solution at infinity vanishes. 
This leads to initial conditions, /s(ri) = 1, fg{ri) = —m? /A. 

Below we present in Table 1 the first 10 states using this shooting method. The 
spin-2 equation is equivalent to the dilaton equation solved in Refs. and 0, so the 
excellent agreement with earlier values validates our method. In all of the equation we use 
a standard Mathematica routine with boundaries taken to be x = — 1 = e and 1/x = e 
reducing e gradually to e = 10~^. Note that since all our eigenfunctions must be even in r 
with nodes spacing in x = — 1 of O(m^), the variable 1/x is a natural way to measure 
the distance to the boundary at infinity. For both boundaries, the values of e was varied 

9 



Glueball by Shooting Technique 



level 


jPC _ 0++ 


JP(^ = 1-+ 


jPC _ 2++ 


n= 


5.4573 


18.676 


11.588 


n= 1 


30.442 


47.495 


34.527 


n= 2 


65.123 


87.722 


68.975 


n= 3 


111.14 


139.42 


114.91 


n= 4 


168.60 


203.99 


172.33 


n= 5 


237.53 


277.24 


241.24 


n= 6 


317.93 


363.38 


321.63 


n= 7 


409.82 


461.00 


413.50 


n= 8 


513.18 


570.11 


516.86 


n= 9 


628.01 


690.70 


631.71 



Table 1: Radial Excitations for Glueballs associated with the Gravitation tensor field. 

to demonstrate that they were near enough to r = 1, and oo so as not to substantially 
effect the answer. 



4. Main Features of Glueball Spectrum 

There are two intriguing aspects of our numerical results for the glueball spectrum. 
The first concerns the pattern of low- mass states, e.g., for the respective ground states, 
one has: m| q < m|n q < rriyQ. The second interesting feature is the asymptotic growth for 
the masses of their radial excitations. In order to provide further insights concerning these 
features, we turn next to approximation schemes appropriate for discussing: (1) ground 
state masses, and (2) high mass excitations. We end this section by providing a heuristic 
picture for understanding these main features. 

4.1. Low-Mass Spectrum 

Eq. ( p.l8|) for tensor modes can be expressed in a standard "Sturm-Liouville" form 

(d d \ Tn^ 

- —x{x + l){x^2)— +w{x)](l){x) = —(t){x) , (4.1) 
dx dx / A 

where we have found it more convenient to use x = — 1 as the new variable. For 



tensor modes, the potential term is in fact absent, i.e., wt{x) = 0. Eq. (p.20| ) for vector 
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modes can also be brought into this form, with = qv{r)/y/r^ — 1 and wv{x) = 

l/{x{x+l){x + 2)). As a Sturm-Liouville system finding various excitations can be carried 
out by applying the minimum principle. 

Using the fact that the ground state energy in such a system corresponds to the 
absolute minimum of its energy functional, we can prove rigorously an inequality 

™t,o < ™F,0: (4-2) 

between the masses for the tensor and the vector ground states. This inequality follows 
immediately from the fact that the potential function for vector modes, wv{x), is strictly 
positive over the range < a; < oo. (See Appendix B.) 

We next apply a variational approach to obtain upper bounds for the respective 
ground state masses. Simple trial wave functions satisfying the respective tensor and 
vector boundary conditions at a; = and at a; — > oo are ^t{x) = (a; + and 
^v{x) = {x + l)~^^/x(x~+2). One readily obtains the following simple upper bounds 
for the ground state masses: 

m|_o < 12, and m^^o < 20. (4.3) 

These values compare very well with the precise numerical answers: m^Q = 11.588 and 
my Q = 18.676. One can easily improve these already impressive bounds by using a better 
variational ansatz. For instance, using (j)T{x) = (a; + 1/a; + X)4>t{x) and 0y(a;) = (x + 
1/x + A)0y (a;), with A as variational parameter, we obtain 

m| o < 11-588, for A = 2.0857 , (4.4) 

m^o < 18.687, for A = 1.1729 . (4.5) 

The situation for scalar modes is slightly more involved. Nevertheless, an analogous 
analysis can also be carried out, as explained in Appendix B. For scalar wave- functions, 
05(a;) = (r^ — l)qs{r)/r, the simplest choice for ground-state trial function consistent with 
boundary conditions at a; = and a; — > oo is 0s (x) = constant, which leads to an upper 
bound for the scalar ground-state mass: 

m|^o < 5.5213. (4.6) 

Observe immediately that this upper bound is already very close to the exact numerical 
result: rrig q = 5.4573. An improved variational bound can obviously be obtained, as was 
done for both the tensor and the vector ground states. It is also important to note that 
this simple variational upper bound for the scalar ground state is already much lower than 
the mass for the tensor ground state. 
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4-2. WKB Estimates For High-Mass Excitations 

We next focus on high-mass excitations by carrying out a WKB analysis. We begin 
by first converting our Sturm-Liouville type equations into "radial" type Schroedinger 
equations, 

{-h^-^+Vix;m^))tPix)=EtP{x), (4.7) 

where E 0~ . In keeping with the standard practice of WKB methods we have also 
introduced a formal expansion parameter h, which at the end of the calculation we set 
back to h = 1. The relations between the new set of wave functions, tl^vix), i(^t{x) and 
"ipsix), and the previous set, (pvi^x), (pvi^) and (/>s(a:), are given in Appendix C. Note that, 
in this approach, eigenvalue parameter m? appears as a parameter in an effective potential, 
and it controls the "strength" of "central force" necessary in forming "zero-energy" bound 
states. For all three cases, when m? becomes large, it can be shown that each potential 
approaches a smooth large-m limit: 

V{x-m^) mVo(x) +0(1), 

1 (4-8) 
4x(x + l)(x + 2)' 

By dividing Eq. ( ^4.7| ) by m^, we observe that, in the large m-limit, % and m appear 
naturally in a simple ratio, %/m. Therefore, a large-m limit is formally the same as 
a"smaU"-^ limit. (Strictly speaking, this is not the case for our scalar equation.) Or less 
formally stated, as we observe in our numerical solutions, the mass eigenvalues increase 
monotonically with the number of zeros n in the wave function so that the large eigenvalues 
expansion is also a short distance limit, as expected. 

For bound-state problems, the WKB consistency condition for eigenvalues m = nin 
for n = 0, 1, 2, , is given as an expansion in 

{n + i) TT^ = /o(m) + him) + h{m) + , (4.9) 

where the coefficients hkim) are expressed as an integral between two classical turning 
points over integrand expressed in term of the potential. Since the leading behavior of the 
integrals are l2k{m) ~ m}~'^^ at large m, we obtain an asymptotic expansion in 1/m at 
Tl = 1. Systematically expanding in 1/m, the WKB condition takes the form: 

(n + i) 7r = som + si + — + ■••, (4.10) 
2 m 
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where coefficients Si are easily calculable. 

We have shown in Appendix D that, in order to obtain the first two coefficients, sq 
and si, it is sufficient to consider only the leading order WKB contribution. 



1 + 4e 

/o(m) = J dx\l~V{x;m?) -^-^^ (4.11) 

where Xi(m) and Xfi{m) are two classical turning points with e 0+. The shift from 
V{x;m'^) to y (x; m^) + first discussed by K. E. Langer [^, is commonly known as 
the "Langer correction". It incorporates the correction at the singular boundary when 
there is a centrifugal term, l{l + l)/x^, in a radial wave equation. Furthermore we have 
added an e term, l/4a;^ — (1 + 4e)/4x^, as a simple procedure for properly implementing 
the threshold behavior for i(j{x) at x = 0. (See Appendix D.) 

^From Eq. ( [4.11J ), one notes that Io{m) can be a non-trivial function of m through its 
dependence in V{x; m?) and on the two classical turning points, ^^(m) and xji{m). Since 
V ~ for m large, it follows that /o(?t^) = m(/o,o + "^~^ -^0,1 + -^0,2 + ■ ■ ■ ) at large m. 
Consequently sq = Iq^q, and si = Iq^i- (To obtain S2, one will also need to calculate l2-) 
Truncating Eq.( [l.lCiD to this order and solving for m^, one obtains the following expansion 



for the mass spectrum, |]n 



ml/n'^ ~ + 5 n + 7, (4.12) 

where /U^, d, and 7 are related to sq, si and S2 hy fx = tc/sq ,5 = 1 — 2si/7r, and 
7 = 574 - 2soS2/7r2. 

We have carried out this WKB analysis for all three cases, and have obtained the fol- 
lowing interesting results. First, we have found that the leading coefficient, /U^, is universal, 
i.e., 

2 1 11 

^-' = (-)y^ rfxv/^M^=-i?(-,-), (4.13) 

is the same for tensor, vector, and scalar modes. Therefore, it is also the same as that 
for the dilaton modes, with /i^ ^ 5.742 . Second, for the term linear in n, simple integer 
coefficients are obtained for both the tensor and vector modes: 



Ot = o Oy = "i- 

However, for scalar modes, a more intricate calculation is needed. We find that 



(4.14) 



= ^(1- ^^^^3^) + ^ + 1-3.018. (4.15) 
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where Ur is the positive root of 9^^ + 3u^ — 1 = 0, — 0.392 . Putting these together, we 
have 

m| ^5.742 (n^ + 3n + 7t), 

~5.742 (n^ + 4n + 7v/), (4.16) 
ml ~5.742 (n^ + 3.018n + 75). 

In Figure 1, we show a comparison on a log- log plot of our WKB prediction to our direct 
numerical results up to n = 9 for all three distinct modes. Interestingly, except for the 
lowest scalar mass, our WKB results agree with the actual data extremely well, with 
constants, 77- — 2.011, 7^ ~ 3.283, and 75 ^ 1.25 fitted to the asymptotic form at large 
n. Higher order WKB results and other related work will be reported elsewhere. 
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Fig. 1: Comparison of WKB approximation with the levels computed in Table 1. 

In spite of this impressive agreement between the low order WKB calculations and 
the exact numerical results, a cautionary note should be added. As pointed out in Ref 0, 
the question of convergence of 1/m-expansion for the logarithm of wave functions, (pvi^), 
(j)v{x) and ^s(x), remains to be addressed. 
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4-3. Heuristic Picture: 

It is possible to provide a qualitative understanding for both the low-mass pattern, 
m| Q < Q < rriyQ, and the asymptotic high- mass excitations by making use of a quan- 
tum mechanical analogy. For m large, Eq. ( [4 .81 ) formally corresponds to a "screened" 
attractive Coulomb potential, with a "bare" charge Q at the core, 

2 

HI 

Q = -^ + 0{l). (4.17) 

More generally, we have shown in Appendix C that, for all three cases, each potential, 
V{x; m^), behaves near a; ~ as 

y(x;m2) ~ — + ^ + 0(1). (4.18) 

This naturally suggests that we identify = Z(/ + 1) as angular momentum squared and 
Q as the unscreened bare charge at the core. We find that ly = 0, and It = Is = 

Our boundary condition for i/j at x = corresponds to choosing il;{x) ~ x^~^^ over 
~ x~K For Re i > — |, this would correspond to always selecting the "regular" solutions. 
In order to avoid the degeneracy ambiguity for the case where = = — i, we shall in 
what follows always assume that It = Is = + ^ and consider the limit e ^ 0+. 

When we include the "Langer correction" , V V + + e){^), the angular momen- 
tum factor, = 1(1 + 1), becomes = (/ + ^)'^ + e. For a; — 0, we find that 

~2 1 m'-6 

W-^^ Qv- 

L| = e, Qt = -^^^, • (4-19) 

~2 ^ — 2 + 16m~2 
^5 = Qs = 3 



At large m^, ( [4.19|) leads to (|4.17|) , in agreement with Eq. (|4.8|) . More generally, since 
V{x; m?) m?VQ{x) for all three cases, this explains our finding that the coefficient /x^ is 
universal. The next order correction, 5, will depend on how V{x;m?) deviates away from 
rin?Vo{x), both at x ~ and at x ~ oo. Indeed, as we show in Appendix D that our results 
for 5t, 5v and 5s follow directly after we have identified the leading asymptotic behavior 
of y at X = and at x = cxo. 

Turning next to the low-mass spectrum. In this analog picture, the vector ground state 
is more massive since it corresponds to a higher angular momentum state when compared 
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to the tensor ground state, Ly > L|, and therefore a stronger binding, (larger m^), is 
required in order to form a zero-energy bound state. In contrast, the tensor modes and 
the scalar modes have the same angular momenta, = L| = e > 0. We find that the 
scalar "potential-well" is "deeper" than that for the tensor modes, especially when is 
small. Indeed, whereas the Coulomb core for tensor modes only become attractive when 
> 12, it is always attractive for the scalar modes, with minimum attraction at = 4. 
Therefore, less binding is required in forming a zero-energy scalar bound state, and it is 
consistent with the fact that the scalar ground state is lighter than the tensor mode. 

5. Discussion 

We have computed the discrete spectrum of metric perturbations in an AdS^ black 
hole background. We would like to end by briefly discussing our results and making a few 
cautionary remarks on the difficulties of the comparison of this strong coupling spectrum 
with glueball masses as computed in lattice gauge theory or classified in bag models. 

As stressed in Ref. [Q, in order to obtain these discrete modes, a non-singular bound- 
ary condition must be applied at the horizon of the black hole and requiring normalizibility 
of wave-functions alone is insufficient for such a purpose. A deeper understanding on the 
physics of the boundary condition used would be highly desirable. Our discussion has 
also been restricted to modes that are singlets under the S0(6) symmetry and the "finite 
temperature" Kluza Klein U(l) gauge group. As a guide to finding states that might 
survive the dimensionally reduced limit, concentrating on states that are non-singlets in 
SO{6) X U{1) appears reasonable. 

Earlier comparisons for QCD^ have noted an apparent agreement between the ratio, 
m{0~~)/m{0~^~^) = 1.45 ± 0.03, found in weak coupling lattice simulations and the ratio, 
m(0~~)/m(0"'""'") ~ 1.50, found at strong coupling from supergravity, using the lowest 
state of the pseudoscalar mode {B^^) and the lowest state of the dilaton ((^). Our new 
lower mass O"*""*" state will somewhat upset this quantitative correspondence. With a much 
smaller value for the scalar mass, this supergravity mass-ratio becomes 3.19 which is too 
large by a factor of two. 

However, the earlier agreement was inexplicably good in any case. From the point of 
view of a 3-d Yang-Mills theory at weak coupling the lowest O"*""*" state should couple to 
both the operator 5" = Tr[F^j^(x)] and the helicity-0 component of the energy- momentum 
tensor T^^, = Tr[F'^^{x)\ — ^rj^^S, which are associated with the dilaton and the scalar 
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component of the metric fluctuations respectively. Indeed each of these operators have 
contributions from the r-component of Af^ along the compact direction, which from the 
view point of 3-d Yang-Mills theory is an adjoint Higgs field. The state associated with 
the Higgs mode must decouple at weak coupling if there is a pure Yang-Mills fixed point. 

A second feature that has previously been noted is the similarity in the level spacing 
of the excitation spectrum of supergravity states and the spectrum of the excited glueball 
spectrum in lattice simulations. In our view, this is even more problematic. For one thing, 
we have so far ignored Kluza-Klein excitations in the compact direction that have a 
level spacing given by Mkk = = 26, where R is the size of the compactified circle 
and the parameter b so defined has units of mass. This spacing is of the same order as the 
mass gap for the lowest 0"''"'" state: m(0'^'^) ~ VS. 457 b ~ 2.3 b. This suggests that the 
proper identification of the excitation spectrum of the supergravity modes with states in 
Yang-Mills theory is non-trivial. In quarkless QCD, all mass splittings are controlled by a 
single mass scale, Aqcd in 4-d or g^N in 3-d, which one may equally well take as the string 
tension, a. In the strong coupling limit on the lattice, this scale is replaced by the UV 
cut-off (or inverse lattice spacing a~^). In these units both the lattice glueball masses and 
the string tension diverge as log{g^N) or log{aglN), whereas in the supergravity scenario 
at strong coupling the putative 3-d "glueball" masses go to constants while AdS string 
tensions diverges as 

(5.1) 

in leading order in strong coupling. It may be that lattice glueball excitations are transverse 
stringy fluctuations that would not survive in the the supergravity approximation. In that 
case, from a lattice view point, there should be no infinite tower of "glueball" states. The 
supergravity spectrum describes radial (5-th dimensional) excitations that either do not 
survive the weak coupling finite N limit or are mixed with stringy states. It is difficult to 
see what should be the precise comparison between these two sets of excitations since they 
have quite different physical features. We consider the understanding of the correspondence 
of this tower of states with Yang-Mills theory at the boundary to be an important unsolved 
issue. 

In spite of these reservations, we believe that a comparison between the supergravity 
results and lattice as well as bag model results does offer opportunities to uncover phe- 
nomenologically interesting features. For instance, in the mass relations for the lowest 
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states as a function of their quantum numbers, we do see a pattern reminiscent of lattice 
glueball spectra. Indeed the pattern, 

m(0++) < m(0++) = m(2++) < m(l"+) 

is generic to both lattice calculations in 3-d and bag calculations, which are in excellent 
agreement with each other. We find that the classification of states with low-dimension 
color singlet fields in a bag picture for glueballs, i.e., in terms of "valence-glue", follows the 
same pattern as the expected in the AdS/CFT correspondence. A fuller understanding of 
these relationships is certainly worth pursuing. First it is evident from our analysis above 
that all 128 lowest bosonic modes for the supergraviton multiplet that are SO{6) x U{1) 
singlets need to be found. For 3-d Yang-Mills, it may be also more fruitful to make 
comparisons with high temperature lattice QCD4 results, where the constant KK adjoint 
field in the compact r direction still remains. 

Our current exercise needs to be extended to schemes for 4-d QCD such as the fi- 
nite temperature versions of AdS'^ x 5"*. As has been suggested elsewhere, the goal may 
be to find that background metric that has the phenomenologically best strong coupling 
limit. This can then provide an improved framework for efforts to find the appropriate 
formulation of a QCD string and for addressing the question of Pomeron intercept in 
QCD. P In addition, a more thorough analysis of the complete set of spin-parity states 
for the entire bosonic supergravity multiplet and its extension to 4-d Yang-Mills models is 
also worthwhile. Results on these computations will be reported in a future publication. 
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August 1999). While this paper was being finalized, a paper [jl5[ appeared, which has a 
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Appendix A. Field Equation for Spin-0 Perturbation: 

Let us start with the gauge where the nonzero components of the spin-0 perturbation 

are 

/iii=/i22=?i(r)e-"^^^ 

h^, = -2Qi(r)^e-"^^3 + g2(r)e-"^"^ (A.l) 
The Einstein equations (|2.10|) give 

93 = - _ 1)2^(3^4 _i) [2^2(0 - tHt' - l)?2,.(r) + A{r' - l)Qi(r)]e— ^ (a.2) 

and one first order and one second order equation in qi and 

4r(3r4 - + ?,r^{r^ + 1)^2,^ - 8(3r4 + l)gi 

_[6r4(r4 - 1) + ^2 r^(3r"-l) ]g^ ^ q ^^'^^ 

— 1 



3r^(r^ - + r[3(r* + 3)(r^ + 1) - Wr^]Qi,^ 

2 

+ [-12(r'^ + if + mV2(3r^ - 1)]qi - m^r^lm^^ — + 3]q2 = 



(A.4) 



r4 - 1 

We can find q2 from (|A.4|) and substitute it into ( |A.3|) to obtain a third order equation for 

qi- 

r3(r4 - 1)2 [mV^ + 3(r4 - l)]qi,rrr 

+r^{r^ - l)[-m^r^ + 2m^r^{r^ + 5) + 3r2(r^ - lf{3r^ + 5)]qi rr 

(A.5) 

+r[-4m^r'^+mV2(23 - 6r^ - r^) - 3(r'^ - l)2(3r^ + 13)]qi,^ 
+ [-mV^ + m'^r^(r^ + 3) - 2mV2(13 + 2r^ + r^) + 48(r^ - l)^]^! = 

It may appear strange that we have obtained a third order system (which needs three 
constants of integration) instead of a second order system as in the spin-2 and spin-1 cases 
above. The source of this extra freedom is that the gauge ( |2.13 ) there is a one parameter 
freedom left in the present ansatz: 

2:3 + ea(r)e~'^^^ 

, ma(r) _ (A. 6) 

r -^r + e— f--e "^^^ 
rj[r) 
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with 

.2 



a(r)=exp[| y(2 + y)] (A.7) 
This freedom gives a change in qi proportional to 



= ' / (A.8) 

(l + r2)¥ 

Substituting 

qi = (A.9) 

we get a second order equation 

r'^ir^ - l)2[-mV2 - 3(r4 - 

+r(r^ - l)[-2m^r^ - mV(17r* - 5) - 3(r^ - l)(9r^ - 

-[mV*^ + 4mV(3r^ - 1) + 4mV2(5r^ - l)(2r^ - 1) + 3(r'^ - l)^(15r'^ + l)]gi = 

(A.IO) 

A similar third order system of equations arises in the study of polar perturbations 
of the Schwarzschild metric, in the gauge where the metric is constrained to be diagonal. 
One discovers a combination of fields that 'remarkably' satisfy a second order equation, 
the 'Zerilli equation' But we can see that such a reduction is to be expected because 



here again there is a one parameter family of diffeomorphisms that preserves the gauge 
condition. Consider a metric of the form 

ds^ = k{r)dt^ + k-^{r)dr^ + r'^{de'^ + sin^ Od^'^) (A.ll) 

Then if we are looking at perturbations with frequency uj in time t, then we have the 
following diffeomorphisms which leaves the metric within the chosen gauge (off diagonal 
components zero): 

t - uh{r)Pi{e)e''^ 

where for infinitesimal e 

h{r) =erk{rf/^. (A.13) 

{Pi(9) is the associated Legendre polynomial in cos 6*.) We can use this solution to reduce 
the third order system to a second order equation. 
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But for the Schwarzschild metric one can choose a different gauge, where such a 
residual gauge freedom does not exist. Such a gauge was used in [0. We thus look for a 
gauge in our problem which would also have no residual gauge symmetry, and so directly 
yield a second order equation. 

Consider the gauge where the nonzero components of the perturbation are 

hrr = ?2(r)e-"^"^ (A.14) 
Ks = hsr = q^{r)e-^-^ 
The Einstein equations (|2.10|) determine q2 and qi in terms of ^3: 



r5(r^- l)gi,^ + 2r^gi 
- ~ (r4-l)2(3r4-l) ^^-^^^ 

qi = ;;^;:5[;J^^^^ (A-16) 

and give a second order equation for: 

r'^{r^ - l)^(3(r4 - l)+m'^r'^)q:i^rr + r{r^ - l)[3(r^ - l){hr'^ + ?,) + m^r'^{7r^ + h)]q^^r 
+ [9(r'^ - l)^+2mV2(3 + 2r^ + 3r^) + mV(r^ - l)]q^ = 0. 

(A.17) 

This is Eq. ( |2.23| ), where we have re-labeled qs by qs- This second order equation is related 
to the second order equation ( |A.10| ) by the transformation 

Appendix B. Sturm-Liouville Approach and Variational Analysis: 

Eq. ( [4.1| ) is in the standard Sturm-Liouville form , 
/ d d \ 



where r(x) = x{x + l)(x + 2) and a{x) = ^. The set of eigenf unctions, {(^n}, provides a 
complete orthonormal basis: < (j)n\<Pm >= dxa{x)(j)riix)(prnix) = 5m,n- It is well-known 
that solving for this set of eigenfunctions and their corresponding eigenvalues is equivalent 
to finding stationery points of the following "energy functional", 

p^^^ ^ dx[T{x)(j)'{xf + w{x)(j){x)^] 
Jq dxa{x)(/}{x)'^. 
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At each stationery point, (^^5 one has r[(/)„] = m^. Since both r(x) and w{x) are positive, 
T[(p] is bounded from below and is in fact positive. In particular, the square of the mass 
for the lowest state, ttiq = T[(f)Q], is the absolute minimum of r[(/)], and uiq > 0. 

Let us denote r^T/^jf*^] as the energy functional for the tensor and vector modes 
respectively. We shall next apply this formalism to obtain a simple and yet important 
mass inequality: 

Denote the the lowest tensor and vector masses by mv,o and rriTfi respectively. Since each 
corresponds to the absolute minimum of its energy functional, it follows that 

< ml = rT[<P^] < TtH; < m^^o = r^[(/)^] < [</>]. (B.3) 

Alternatively, to^q "^^^ expressed as to^q =< >o +rT[<^o']5 where 



< wv >o = 



dxa{x)(j)Q (xY 



Note that wvix) = l/(x(x + l)(x + 2)) is a positive function for < x < oo; it follows 
that < Wv >o is positive, and 

ml,o > rT[<l>^]. (B.4) 

On the other hand, since q is the absolute minimum of Ft, we also have the following 
inequality, 

Tt[<P^] > Tt[<PI] = ml,. (B.5) 

^From ( p.4| ) and (|B.5|) , the mass inequality between vector and tensor ground states, 
rriyQ > m^Q, Eq. (^), follows. 

It is also worth pointing out that, in re-deriving these eigenvalue equations from the 
energy functional, it is necessary to impose boundary conditions 

T{x)(i){x)(t)'{x) 0, (B.6) 

both for X — s> and x — > oo. This states more clearly what boundary condition is to be used 
at the horizon of the black hole. In particular, this rules out the solution where 0(x) ~ log x 
at X = 0. To be more precise, by examining our differential equations, Eq. (|4.1| ), both 
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at x = and at x oo, Eq. ( |B.6| ) requires that tensor and vector wave- functions must 
behave as 

(B.7) 

(prix) ^ constant, (pvix) ^ x^''^ , x ^ 0. 
We next address the question of scalar modes. Let us first convert Eq. ( |2.18|) into 

Lscpsix) = (- (^^5(x;m2)A) + wsix;m^)'^(l)six) = 0, (B.8) 

where 

x(x + 2) , m^\(3x'^ + 6x + 4) -m^(x + l)] 

Ts{x;m') = j-——^^ , ws{x;m')- 



[3x{x + 2)+m^{x + l)y ' ^ 4(a; + l)[3x(x + 2)+m2(a; + l)]2' 

(B.9) 

Here, (psi^) is related to qsir) in (|2.18| ) by (ps{x) = ^^lH^ Qsif)- Since mP enters into 
the differential equation in a rather non-trivial fashion, we are unable to directly convert 
this into a Sturm-Liouville problem with as eigenvalues. Let us consider instead the 
following eigenvalue problem: 

Ls<Pnix) = Enas{x)(j)n{x), (B.IO) 

where Eq < Ei < E2 < ■ ■ ■ , and asix) is at this point unspecified. This in turn can be 
obtained via a Sturm-Liouville type variational approach, i.e., finding stationery points of 
the following energy functional: 

p^^^j ^ /o°° dx{Ts{x;m'^)\(P{xy\'^ + wsjx; m^)\(f){x)\'^) 

dxas{x)\(p{x)\'^. 

Note that, here, is to be treated as a parameter. To have a properly defined stationery 
conditions, we require ts{x; m'^)(f){x)(j){xy = 0, both at a; = and x — 00. This, in turn, 
fixes our boundary conditions. A careful examination indicates that the desired boundary 
conditions for (psix) are 

(ps{x) — > constant, (B-12) 

for both X — >■ 00 and x — 0. 

Since ws{x;m?) is no longer positive, neither is r5[(/)]. Nevertheless, Vs[(p] is still 
bounded from below since ws{x; ni?) is. As such, we identify E'^ijm?) = rs[(pn], in increas- 
ing order, at each local minimum, (p^. In particular, Eo{m'^) is the absolute minimum of 
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r^, although it no longer is positive in general. This "generalized Sturm-Liouville prob- 
lem" reduces to our scalar problem if a zero-energy solution exists. To be more precise, one 
needs to answer the following question: "For what values of m^, j = 0, 1, 2, ■ ■ -, will there 
be a zero-energy solution where Ej{m'j) = ?" As a minimum problem, we are searching 
for uiq for which the absolute minimum of Ts{(p) vanishes, i.e., minr5[^; mg] = 0. Turning 
this into a variational problem, we require that, for > rriQ, 

Q _ ^ ir dx{Ts{x-m'^)\(l){xy\'^ + wsjx; m^)\(j){x)\'^) 

dxas{x)\(p{x)\'^ 

for suitably chosen trial wave- functions. For each choice of (p{x), the value of at which 
the integral on the right vanishes provides an upper bound for uiq. In arriving at our 



bound for the scalar ground-state mass, Eq. (|4.6|) , the trial function (j)s{x) = constant, 
was used. 



Appendix C. Radial Schroedinger Representation and Effective Potentials: 



In order to bring our ODE, Eq. into the radial Schroedinger form, Eq. (| 

we perform the following transformation: 

i^rix) ={x{x + l){x + 2))^ (prix), 

ijv{x) ={x{x + l){x + 2))^(j)vix), ^C.l) 
i^s{x) =[ 3 )^(j)s{x). 

The resulting effective potentials for all three cases are: 

i 2t^ / ^ , 3(x2+2x + 2) 1 



Ax{x + lY{x + 2) x2(x + l)2(x + 2; 



2 ■ 



Vy{x-m') =mVo(x) + 3(^^ + 2^ + 2) ^C.2) 

4x(x -I- l)^(x -I- 2) 

T. / 2^ 2T./N m2(3x2+6x + 2) 

Vs{x- m ) =m V,{x) + + + + 1) + 3,(, + 2)] + 

where a rather involved expression has to be added for the scalar potential, 

U , 1X2 , 1N2 , . / , o^ ^4 I 4m^x(x+l)(x+2)(3(x+l)^+7) 

_ A{x + iy + iy + 4x{x + 2) ./ rn + (((^+1)2+1)2+4^(^+2)) \ , . 

^ ^ 4x2(x-f2)2 ^\ [m2(x + l) -F3x(x + 2)]2 ^^-^^ 
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Let us briefly comment on these effective potentials, V{x;m'^) = m^) + (l/4+e)(l/x^), 
in various limits, 

(a) — > oo: For all three effective potentials, 

V{x-m^) =m'^Vo{x) +0(1). (C.4) 

(b) X — > 0: In going from V to V, for x ~ 0, we have V{x;m'^) ~ I? jx^ + Qjx + 0(1), 
where = (/ + + e. Those and Q values given in Eq. ( [4.19|) follow from 

M-;m^)=^-!^+0(x«), 

yW-;-^)~-^+o(x«), (C.5) 



Vs[x- m^) =4 - ^"^ ; ' + 0(x°). 
x^ 8x 



(^2 _ 2 + 



(c) X ^ oo: 



VT(x;m2)=^-^^ + 0(x-4), 

X^ X'^ 

n^(x; m^) =1 - + 0(x-4), (C.6) 

From these, one concludes that each "Coulomb" bare charge, Q, is always screened at large 
distance; an effective potential vanishing as Xjx^ at large distance is overall "neutral". 



Appendix D. WKB Treatment: 

In a standard WKB treatment for a one-dimensional quantum mechanics problem, 
one relies on having a small parameter, "?i", which allows a systematically expansion 
for the logarithm of the wave-function, (^(x) ~ exp[±|- dx^jE — V{x) -\- ■■■]. For a 
problem involving a spherically symmetric potential, V{x), one has a one-dimensional 
radial Schroedinger equation involving an effective potential, e.g., Eq. ([4.7|), where V{x) 
Veff{x) = V{x) + l{l + l)fb^/x'^. In solving the radial equation, normalizability requires 
that the wave- function if; ~ x'"*"^ as x — O"*", which also corresponds to our "non-singular" 
boundary condition at the black hole horizon. Note that this exact threshold behavior, 
ip{x) ~ g('+i)iogx^ would normally appear only in the next to leading order in a WKB 
expansion. 
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It is indeed possible to arrange so that the threshold behavior is properly implemented 
in the leading order by working with a "rapidity" variable, y = logx, where y is now 
unrestricted, — oo < y < oo, and the radial Schroedinger Equation becomes ( — cP / dy"^ + 
U{y))'4'{y) = E'ip{y), where i>{y) — e~^'i(j{x), E = Ee'^^, and 

U{y) ^ e'yV.ff{x) ^ x^V.ffix) + ^) = x\V{x) + + 

Note that, for U{y) and Veff{x), the combination /(/ + 1)^^ has become (/ + ^)^^^, which 
correlates precisely with the appropriate threshold behavior for ip ~ x^~^^ . For an attractive 
potential V{x) which is less singular than at the origin, this new effective potential U{y) 
will have an attractive well, with U{y) (/ + ^)^^^ when y —oo, whose bound state 
energies can be calculated via a standard WKB approximation. For instance, the leading 



order WKB contribution is {n + ^)n?i = dyy E ~ U{y), where n = 0, 1, 2, ■ ■ -, e^^ = xl 
and e^^ = x^. 

Returning to our problem; with E 0~ , and 1/m playing the role of h, we arrive at 

1 rVRim) 

/o(m) = - / dyy/-U{y;m). (D.l) 

^ JyL{m) 

By changing variable from y back to x, Eq. ( |D.1D can also be written as Eq.( [4.11] ) where 
where Veff{x;m) enters. Let us turn next to the expansion of Io{m) in Io{m) ~ 
-^0,0 + ^-^0,1 + 7)^-^0,2 H — ■■ Before proceeding, we separate the integral in Eq. ( p.l| ) , into 
two, 

loim) = him) + Inim) = — dy^-U{y-m) + — / dy^-U{y-m). (D.2) 

This separation at y = is of course arbitrary; our final result will be independent of 
this choice. It is convenient to separate U{y;m) into two pieces, U{y;m) = m^(C/'o(2/) + 
m~^C/i(y; m'^)), which will be used in Eq. ( p.2|) . To expand in 1/m, one needs to know 
how Uq and Ui behave as y ±oo. 

Let us begin by treating the case of tensor modes with reasonable amount of details. 
Consider /^(m) and let us examine the lower integration limit, With both e ~ 0"*" 

and 1/m 0, the left-turning point is determined by the asymptotic limit of U{ii;m?) 
as y ^ -oo, U{y;m'^) ^ e - - 3/2)6?'. That is, ~ ^J^- As e ^ 0+, for m 

large, it follows that the left-turning point moves to to minus infinity, yL — oo. More 
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explicitly, this is a consequence of the fact that = e ^ 0. We further note that, with 
both t/o(y) and Ui{y) vanish as O(e^) as y ^ — oo, the ratio Ui{y)/UQ{y) remains finite; it 
is justified to expand the integrand in and one obtains 

lL{m)= f dy^^U^) + Q{\) = [ dx./^Vo{x)+0{^). (D.3) 

J-oo ^ Jo ™ 

We turn next to lR{rn), which can be more conveniently written as 



= / -7=in^+y^ _y=inw^){\/i + ;^-^-i} (D.4) 

The upper limit is determined by V{xYi) = m^(Vo(xi?) + ;;^Vi(xij)) = 0. Since 
Vo{x) — > —l/4x^ and Vi{x) — it follows that XR{m) ~ m? /A. The first integral 
can be expressed as a difference: j'^ dx\J —Vq{x) — dx^J—VQ{x), where, for 1/m small, 
dx\J —Vo{x) ~ ^R^"^ — 2/m. For the second integral, the main contribution now 
comes from region where x = O(m^). Let us change variable from a; to w = x/x^ ~ Ax/m?, 
and, in the limit of large m, this expression becomes ^ ;jjy2 (Vl — u — l) + O(^) = 
(^)(1 — f ) + 0(;^). Putting things together, we thus have, 

/■oo -1 

lR{m)= dx^/-Vo{x)--+0{^). (D.5) 
Combining this with that for 7^, we have, for tensor modes. 



I(m) 



poo ^ 1111 1 

/ dx,/-Vo{x) -- + 0( — ) = -B{-, -)-- + 0( — ). (D.6) 
Jq m 4 4 2 m 



This leads to Eq. (|4.13|) for /i, and 5t = 3 for tensor modes. 



Turning next to the vector modes. Again, separating /(m) into two parts, Eq. ( p.2|) , 
an identical analysis for lR{m) can be carried out as was done for the tensor modes, leading 
to the expansion, Eq. ( p.5| ). However, since U{ii) 1/4 as y — oo, (corresponding to 
the situation where Ly = 1/4), Eq. ( p.3|) no longer holds for the vector modes. One finds, 
instead, 

lL{m)= [ dx^/^Vo(x)-^ + 0{^). (D.7) 
Jq 2m 

Combining this with that for Ir, we have, for vector modes. 



/(m)^/ rfxv/^M^- — +0( — ). (D.^ 
Jq 2m 
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This leads to Eq. ( [4.13|) for n, and Sy = 4: for vector modes. 

We complete this analysis by treating the scalar modes. Since L| = Lj, = 0, one 
finds that Eq. (|D.3|) holds for scalars modes also. As for Iji{m), we need to work with 
Eq. ( p.4| ) with V^x^rn?) = l/(4x^) + y(x;m^), where V for scalar modes is given by the 
rather involved expression in ( P.2| ) and ( P.3|) . Let us consider the second integral first, 
where the main contribution now comes from region where x ~ rm? . Changing variable 
from x to w = x/m?, and, in the limit of large m, this expression becomes 



1 r« du .Vl-3M2_9^i3 1 

where ur is the positive root of the 1 — "iv? — = 0, ur ~ .3915 . The first integral 
can be expressed as a difference: /^°° dx^J —Vq{x) — dx^J—VQ{x), where, for 1/m small, 
dx^f—Voix) ~ l/{m^/uj^. Combining these with that for 1^, we have 

/(m)~ / dx^J-Vo{x) + - + Q{^), 
Jq m 



where si ~ -(2.018) (7r/2), thus leading to 5s = 3.018 , Eq. (|iT5D . 
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